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IDEAL  ELASTIC,  ANELASTIC,  AND  VISCOELASTIC 
FLOW  IN  A  METALLIC  GLASS 

A.I.  Taub  and  F.  Spaepen* 


INTRODUCTION 

When  a  material  is  subjected  to  an  applied 
stress,  its  strain  response  generally  consists  of 
several  components,  each  differing  by  its  depen¬ 
dence  on  time,  by  the  degree  to  which  it  is 
recovered  upon  removal  of  the  stress,  and  by  the 
linearity  of  the  response. u>  Glassy  metallic  alloy 
systems  have  been  reported  to  exhibit  at  least  four 
of  these  strain  components: 

1.  Ideal  elasticity  (recoverable,  instantaneous, 
linear  stress-strain) (2  J) 

2.  Anelasticity  (recoverable,  time-dependent, 
linear  stress-strain)  (4~8’ 

3.  Viscoelasticity  (permanent,  time-dependent, 
linear  stress-strain  rate)'5-6-910’ 

4.  Instantaneous  plasticity  (permanent,  instan¬ 
taneous,  nonlinear  stress-strain)"1’ 

In  this  study,  an  amorphous  Pd-based  alloy  was 
tested  in  the  homogeneous  flow  regime.  This  elim¬ 
inates  the  strain  contribution  of  the  instantaneous 
plastic  flow  mechanism. <1J  U)  For  each  of  the  three 
remaining  strain  components,  the  time  dependence, 
recoverability,  and  linearity  are  examined  in  detail. 
The  effect  of  structural  relaxation  on  each  type  of 
strain  is  also  discussed. 

RESULTS 

The  tensile  creep  flow  of  amorphous  PdS2Siu 
wires  was  measured  by  the  technique  described  pre¬ 
viously. <91  Both  as-quenched  and  preannealed  speci¬ 
mens  were  tested.  The  preanneals  served  to  stabi¬ 
lize  the  structure  of  the  specimens,  thus  permitting 
isoconflgurational  testing.  For  preannealed  speci¬ 
mens  that  were  to  undergo  stress  reductions  during 
the  creep  test,  load  cycling  was  conducted  during 
the  anneal  in  order  to  eliminate  any  structural  re¬ 
laxation  that  might  be  associated  with  load  removal 
during  the  test."4’ 

Figure  la  shows  a  typical  strain  vs  time  plot  ob¬ 
tained  for  a  preannealed  specimen  after  a  stress  in¬ 
crease.  The  preanneal  in  this  case  was  32S  hours  at 
$00  K.  The  instantaneous  elongation  yE  is  charac¬ 
terized  by  ideal  elasticity.  The  steady-state  strain 
y p,  established  after  approximately  SO  hours,  is 
viscoelastic.  The  anelastic  strain  contribution  y A  is 
the  transient  flow  that  occurs  before  the  steady- 
state  condition  is  established. 

'Hinird  University 
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Figure  1.  Strain  response  to  equivalent  shear  stress 
Increase  from  105  to  156  MPa.  Sample 
preannealed  at  testing  temperature  of  500  K 
for  343  hours.  Symbols  x  and  a  are  actual 
data.  Solid  lines  are  least-square  curve  tits 
using  the  five  exponential  decay  formulation 
described  in  the  test. 

(a)  Total  equivalent  shear  strain  (y  -  s? 
(•uniaxial  tensile  strain) 

(b)  Equivalent  anelastic  shear  strain  y ,, 
obtained  by  subtracting  the  ideal  elas¬ 
tic  and  viscoelastic  strains. 

Ideal  Elasticity 

The  eLstic  response  of  specimens  preannealed 
at  500  K  for  325  hours  was  measured  at  293,  424 
and  500  K.  In  all  -cases,  the  response  was  instan¬ 
taneous,  completely  reversible,  and  linear. 


The  tensile  elastic  modulus  (i.e..  Young's 
modulus)  at  293  K  was  determined  to  be  84 
±8  GPa.  This  is  in  good  agreement  with  Davis’00 
value  for  PdgoSijo  of  88  GPa. 

Using  the  values  of  the  moduli  determined  at  all 
three  temperatures  and  assuming  a  linear  tempera¬ 
ture  dependence,  we  found  that  d( In  £/  dT  *= 
(—3.2  ±0.6)  x  10-4C-1.  This  is  to  be  compared 
with  the  P'-ry  and  Pritchet05)  dynamic  measure¬ 
ment  of  —2.9  x  10-4C-1. 

Structural  changes  in  the  amorphous  state  have 
been  shown  to  affect  the  elastic  stiffness  of  metallic 
glasses.0’  Increases  in  modulus  approaching  10%, 
relative  to  the  as-cast  condition,  have  been  reported 
for  many  systems  annealed  near  the  glass  transition 
temperature/3’  In  our  tests,  the  32S-hour  preanneal 
at  300  K  was  sufficient  to  stabilize  the  structure, 
and  no  changes  in  modulus  were  observed  during 
the  loading  and  unloading  tests. 

Anelastlcity 

Figure  lb  shows  the  anelastic  component  of  the 
total  strain  of  Figure  la,  obtained  by  subtracting  the 
ideal  elastic  and  viscoelastic  contributions.  The 
time-dependent,  transient  nature  of  the  flow  is  evi¬ 
dent.  Additional  anelastic  flow  is  not  resolvable 
after  approximaieiy  two  days. 

Complete  recoverability  is  a  requirement  for 
true  anelastic  behavior.  In  Figure  2,  the  results  of 
a  test  of  this  condition  are  examined  by  observing 
the  response  to  a  stress  cycle  from  35  to  71  to 
33  MPa.  The  entire  response  is  plotted  in 
Figure  2a.  The  elastic  response  is  seen  to  be  com¬ 
pletely  and  instantaneously  recovered  on  stress 
reduction  (y£—  —yE').  In  Figure  2b,  only  the  an¬ 
elastic  contribution  is  shown.  Complete  recovera¬ 
bility  of  this  strain  response  is  verified 

y a  “  -?  A 

The  linearity  condition  for  true  anelastic  flow  re¬ 
quires  that  the  total  anelastic  strain,  after  the  com¬ 
plete  decay  of  the  transient,  be  directly  proportional 
to  the  magnitude  of  the  stress  change.  To  check 
this  requirement,  a  creep  sample  was  subjected  to 
sequential  stress  increases  at  500  K,  and  the  anelas¬ 
tic  strain  was  determined  in  the  manner  illustrated 
in  Figure  1.  In  Table  1,  the  observed  anelastic 
strain  is  listed  for  the  stress  increments,  and  a 
linear  relation  is  verified,  y A/br  =  (8.0  ±0.80)  x 
10-4  cm/cm-MPa,  where  r  is  the  equivalent  shear 
stress  t  —  o7>/3  and  <r  is  the  uniaxial  tensile  stress. 

A  corollary  of  the  linearity  postulate  is  the  su¬ 
perposition  (additivity)  of  the  responses/”  In 
another  load  cycle  experiment,  this  principle  was 
checked.  The  stress  was  raised  in  two  steps  from 


(b) 

Figure  2.  Strain  response  to  stress  cycle  from  3S  to  71 
to  35  MPa.  Sample  preannealed  at  testing 
temperature  of  500  K  for  575  hours. 

(a)  Total  strain 

(b)  Anelastic  strain. 


Table  1 

TOTAL  ANELASTIC  STRAIN 
FOR  VARIOUS  STRESS  INCREASES* 


Equivalent  Shear 
Stress  Range 
r  —  <r/>/3 
(MPa) 

Total  Anelastic 
Equivalent  Shear  Strain 
yA  “  >/3  tA 
(lO-4  cm/cm) 

(10-*  MPa- 

54  to  105 

4.0  ±0.3 

7.8  ±0  6 

105  to  156 

4.7  ±0.3 

9.0  ±0  6 

156  to  207 

4.2  ±0.3 

8.1  ±0.6 

207  to  259 

4.2  ±0.3 

8.1  ±0.6 

259  to  310 

3.6  ±0.7 

7.2  ±1.2 

310  to  363 

3.8  ±0.3 

7.2  ±0.6 

363  to  383 

1.7  ±0.3 

8.4  ±1.8 

’Specimen  preannealed  ai  testing  temperature 
of  $00  K  for  343  hours. 
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35  to  71  to  106  MPa.  The  associated  anelastic 
strains  were  (3.5  ±0.2)  x  10“4  cm/cm  and  (3.8 
±0.3)  x  10-4  cm/cm.  The  stress  was  then  reduced 
in  one  step  to  35  MPa,  with  a  corresponding  anelas¬ 
tic  recovery  of  (7.0  ±0.3)  x  10-4  cm/cm.  The  su¬ 
perposition  of  the  responses  was  thus  verified. 

Although  several  investigators  have  examined 
the  effect  of  structural  changes  on  the  anelastic  re¬ 
laxation  spectrum  of  metallic  glasses, (2-4)  no  sys¬ 
tematic  investigation  of  the  effect  of  structural 
change  on  the  anelastic  creep  response  of  these 
glasses  has  been  reported.  We  observed  no  change 
in  the  anelastic  response  of  the  preannealed  speci¬ 
mens  during  the  testing,  indicating  that  the  325- 
hour  preanneal  at  the  testing  temperature 
sufficiently  stabilized  the  structure.  It  is  possible 
that  the  anelastic  creep  response  changes  during  the 
early  stages  of  annealing,  but  we  have  not  made  ob¬ 
servations  in  this  regime. 

Viscoelasticity 

The  viscoelastic  flow  exhibited  by  well-annealed 
samples  (i.e.,  no  structural  relaxation  occurring 
during  the  test)  is  characterized  by  a  constant  strain 
rate.  Referring  to  Figure  1,  the  viscoelastic  com¬ 
ponent  is  shown  to  contribute  throughout  the  test, 
although  constant  strain  rate  steady-state  flow  is  not 
fully  established  until  after  the  decay  of  the  anelas¬ 
tic  transient. 

Viscoelastic  flow  must,  by  definition,  produce  a 
permanent  strain.  The  test  shown  in  Figure  2a  can 
be  used  to  check  this  requirement.  The  total 
viscous  flow  yp,  up  to  84  hours  (when  the  stress 
was  reduced),  is  (5.0  ±0.3)  x  10~4  cm/cm.  The 
observed  permanent  set  yp'  is  (5.7  ±0.5)  x 
10~4  cm/cm,  verifying  nonrecoverability  of  the 
flow. 

The  linearity  of  the  viscoelastic  stress-strain  rate 
relation  has  been  discussed  elsewhere. <l6>  A  transi¬ 
tion  from  linear  to  nonlinear  behavior  was  shown  to 
be  an  inherent  property  of  metallic  glasses.  For 
PdS3Sit8  specimens  preannealed  and  tested  under 
the  same  condition  as  those  discussed  in  this  study, 
the  limiting  condition  for  linear  flow  was  esta¬ 
blished  at  an  equivalent  shear  stress  r  -  <r/V 3  = 
300  MPa,  where  cr  is  the  uniaxial  tensile  stress. 

The  temperature  dependence  of  the  iso- 
configurational  viscoelastic  strain  'ate  has  been 
measured.  The  results  are  reported  in  terms  of  the 
viscosity  ij  -  r/yp  (see  Figure  3).  To  stabilize  the 
structure,  the  specimens  were  preannealed  under  an 
argon  atmosphere  with  no  applied  stress.  The  an¬ 
nealing  temperatures  and  times  are  indicuted. 
These  preanneals  enabled  us  to  test  the  specimens 


Figure  3.  Isoconfigurational  viscosities  obtained  after 
the  indicated  annealing  treatments. 

in  the  temperature  range  490  to  525  K  without  ob¬ 
serving  additional  structural  relaxation.  The  ab¬ 
sence  of  structural  change  was  verified  in  each  case 
by  returning  to  the  initial  testing  temperature  at  the 
end  of  the  temperature  cycle  and  observing  no 
significant  change  from  the  initial  value  of  the 
viscosity.  The  activation  energies  found  for  the 
isoconfigurational  viscosity  of  these  specimens  is 
200  ±15  kJ/mole,  in  agreement  with  other  mea¬ 
surements  for  specimens  preannealed  at  lowe>  tem¬ 
peratures. (9>17)  Tests  at  different  stresses  also 
showed  no  change  in  the  activation  energy. 

Although  the  activation  energy  for 
isoconfigurational  flow  is  stress  and  thermal  history 
independent,  the  magnitude  of  the  viscoelastic 
strain  rate  is  highly  dependent  on  the  degree  of 
structural  relaxation,  and  therefore  on  thermal  his¬ 
tory.  Previous  investigations*9-17’  have  shown  that 
the  viscosity  can  be  changed  by  many  orders  of 
magnitude,  even  when  annealing  at  temperatures  as 
low  as  Tt- 200.  Figure  4  shows  viscosity-time  his¬ 
tories  for  specimens  tested  under  different  stresses 
at  500  K  from  the  cast  state  (curves  A,  B,  and  C). 
The  linear  increase  of  viscosity  with  time  has  previ¬ 
ously  been  observed. n7-,*)  Note  that  for  the  stress 
range  tested  (39  <  r  <  155  MPa),  the  viscosity 
annealing  kinetics  arc  stress-independent.  Further- 


3 


Figure  4.  The  viscosity  as  a  function  of  annealing 
time  for  specimens  tested  at  the  indicated 
stresses.  Samples  A,  B,  and  C  were  creep 
tested  from  the  as-cast  state.  The  stress  on 
sample  A  was  increased  after  80  hours  of 
testing.  Sample  D  was  preanneaied  as 
described  in  the  text.  The  origin  for 
sample  D  was  shifted  to  t  «■  50  hours. 

more,  there  is  evidence  that  subsequent  annealing 
kinetics  are  not  affected  by  previous  thermal  histo¬ 
ry.  Curve  D  is  for  a  specimen  that  was  sequentially 
preannealed  at  424  K  for  225  hours,  444  K  for 
175  hours,  466  K.  for  152  houis,  and  487  K  for  240 
hours.  The  specimen  was  then  creep  tested  at 
502  K.  During  the  test,  the  viscosity  was  found  to 
increase  linearly  with  time  from  an  initial  value  of 
5  x  1015  Ns/m2.  In  Figure  4,  the  vi  cosity-time 
plot  for  this  test  is  shown  with  the  origin  shifted  to 
coincide  with  the  other  curves.  The  subsequent  an¬ 
nealing  kinetics  are  seen  to  be  the  same  as  those 
for  the  as-cast  specimens,  implying  that  the  kinetics 
of  structural  relaxation  are  independent  of  the  ther¬ 
mal  history. 

For  certain  properties,  such  as  the  Curie  tem¬ 
perature,09'  a  “cross  over  effect”  is  observed  dur¬ 
ing  structural  relaxation,  indicating  a  thermal  histo¬ 
ry  memory.  On  the  other  hand,  the  viscosity  does 
not  appear  to  exhibit  this  property. 

DISCUSSION 

Anelastic  flow  in  many  crystalline  alloys  can  be 
characterized  as  a  single  process  with  the  general 
form  of  an  exponential  decay:  y  A  -  yJO— 

We  attempted  to  fit  the  flow  observed  in  these  tests 
with  this  relation  combined  with  a  viscoelastic  term 
yO)-y2<l-  e~,/i)  +  yFr.  A  typical  least-squares 
fit,  using  y\,  r,  and  yF  as  free  parameters,  is 
shown  in  Figure  5.  Only  the  anelastic  strain  is 
shown,  the  viscoelastic  strain  y  Ft  having  been  sub¬ 
tracted  from  the  total  strain.  As  a  resuit  of  this 


Figure  5.  The  equivalent  anelastic  shear  strain  com¬ 
puted  from  the  data  of  Figure  la,  using  a 
single  exponential  decay  formulation.  The 
symbols  are  the  actual  data  points,  from 
which  the  fitted  viscoelastic  strain  has  been 
subtracted.  The  solid  line  is  the  least- 
squares  fit. 

subtraction,  the  anelastic  data  peaks  at  approximate¬ 
ly  70  hours  and  then  drops.  This  is  physically  un¬ 
realistic,  and  is  only  an  artifact  of  the  fit  value  for 
the  viscoelastic  strain  rate  yp  being  higher  than  the 
actual  value.  In  addition,  the  data  starts  higher 
than  the  fined  curve  and  then  crosses  it  twice  be¬ 
fore  peaking.  Similar  systematic  deviations  were 
observed  in  every  case  examined  in  this  manner. 
This  incompatibility  of  the  data  with  a  single  pro¬ 
cess  exponential  decay  is  an  indication  of  the  ex¬ 
istence  of  an  anelastic  site  spectrum.  The  existence 
of  such  a  spectrum  in  the  more  common  glassy  sys¬ 
tems  is  well  established.'20'  Recent  work  on  metal¬ 
lic  glasses  has  shown  that  anelastic  flow  in  these 
materials  is  also  governed  by  a  spectrum.015’  The 
anelastic  flow  can  then  be  expressed  as 
yA  -  yJA(\-e~'lrj),  where  the  sum  is  replaced 
by  an  integral  for  the  case  of  a  continuous  spec¬ 
trum.  Fitting  experimental  data  by  an  exact 
method  to  such  a  relation  is  difficult,  and  one  must 
resort  to  an  approximation  technique. l2,,  We  chose 
the  finite  spectrum  approach,  employing  a  least- 
squares  fit  to  the  prefactor  of  five  exponential  de¬ 
cays,  with  predetermined  time  constants  spaced  by 
a  factor  of  three  (iy  *  0.2,  0.6,  1.8,  5.0,  and  17 
hours).  The  smallest  time  constant  is  limited  by 
the  data  acquisition  rate,  one  point  every  0.1  hour. 
The  decay  with  the  largest  time  constant,  17  hours, 
reaches  0.95  of  its  final  final  value  after  50  hours. 
From  the  recovery  tests,  it  is  known  that  additional 
anelastic  flow  is  undetectable  after  this  time. 
Therefore,  this  spectrum  spans  the  necessary  time 
constants.  Performing  a  least-squares  fit  to  the  data 


with  this  relation  combined  with  a  viscoelastic  term 
y (<)  -  y^(l-e"'/T-')  +  ypt,  and  using  yj,  and 

yf  as  free  parameters,  we  produced  an  excellent  fit 
to  the  data  (sec  the  solid  lines  in  Figure  1).  Note 
that  the  fit  value  for  the  viscoelastic  strain  rate  yp  is 
reasonable  since  it  does  not  produce  an  artificial 
peak  in  the  anclastic  strain.  This  method  provides 
an  analytic  expression  from  which  the  total  anelas- 
tic  flow  can  be  computed.  However,  fitting  the  data 
with  a  finite  sum  of  exponentials  whose  time  con¬ 
stants  are  arbitrarily  predetermined,  yields  a 
nonunique  solution,  and  therefore  no  physical 
significance  could  be  attributed  to  the  weights  of 
the  resulting  spectrum.0-21* 

CONCLUSIONS 

The  homogeneous  strain  response  of  amorphous 
PdgjSiu  to  an  applied  stress  consists  of  three  com¬ 
ponents:  elastic,  anelastic,  and  viscoelastic  strains. 

The  elastic  response  is  fully  reversible,  instan¬ 
taneous,  and  linear.  The  tensile  elastic  modulus  at 
293  K  was  determined  to  be  84  ±  8  GPa.  The  tem¬ 
perature  dependence  is  d(,lnE)/dT  *  (-3.2  ±0.6) 
x  10-4C-'. 

The  anelastic  response  is  a  transient  that  follows 
stress  increases  and  is  fully  recovered  upon  stress 
reduction.  For  a  Pdg2Silg  specimen,  preannealed 
and  tested  at  500  K,  the  anelastic  transient  contri¬ 
butes  measurably  of  the  flow  for  approximately  50 
hours.  The  transient  cannot  be  described  by  a  sin¬ 
gle  process  exponential  decay.  However,  a  sum  of 
exponential  decays  spanning  a  spectrum  of  time 
constants  from  0.2  to  17  hours  provides  an  excel¬ 
lent  fit  to  the  data. 

The  total  anelastic  strain,  after  the  complete  de¬ 
cay  of  the  anelastic  transient,  varies  linearly  with 
the  magnitude  of  the  stress  change  over  the  entire 
stress  range  tested:  y ,,/At  ■*  8.0  ±0.8  x  10-6 
MPA"*.  The  superposition  principle  was  shown  to 
apply  to  anclastic  flow. 

The  isostructural  viscoelastic  flow  is  nonrecover- 
able  and  characterized  by  a  constant  strain  rate. 
The  activation  energy  for  isoconfigurational  flow 
was  found  to  be  200  ±15  kJ/moIe,  independent  of 
the  annealing  temperature  used  to  stabilize  the 
structure.  When  structural  relaxation  occurs,  the 
viscoelastic  strain  rate  changes  by  many  orders  of 
magnitude.  The  strain  rate  has  been  observed  to 
vary  inversely  with  time.  The  rate  of  decrease  was 
found  to  be  both  stress  and  thermal  history  in¬ 
dependent. 

The  viscoelastic  strain  rate-stress  relation  is 
linear  only  in  the  low  stress  regime,  a  transition  to 
nonlinear  behavior  occurring  at  approximately 


300  MPa.  The  anclastic  strain,  on  the  other  hand, 
remains  linear  over  the  entire  stress  range  tested 
(0  to  383  MPa) 
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